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Abstract 

The longwave phenomenological model is used to make simple and precise calculations of various physical 
quantities such as the vibrational energy density, the vibrational energy, the relative mechanical displacement, and 
the one-dimensional stress tensor of a porous silicon distributed Bragg reflector. From general principles such as 
invariance undertime reversal, invariance under space reflection, and conservation of energy density flux, the 
equivalence of the tunneling times for both transmission and reflection is demonstrated. Here, we study the 
tunneling times of acoustic phonon packets through a distributed Bragg reflector in porous silicon multilayer 
structures, and we report the possibility that a phenomenon called Hartman effect appears in these structures. 
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Background 

Phonons, the quanta of lattice vibrations, manifest them- 
selves practically in all electrical, thermal, and optical 
phenomena in semiconductors and other material sys- 
tems. The reduction of the size of electronic devices 
below the acoustic phonon mean free path creates a new 
situation for phonon propagation and interaction, open- 
ing up an exciting opportunity for engineering phonon 
spectrum in nanostructured materials [1]. Since the early 
work of Narayanamurti et al. [2], important progress has 
lately emerged in the development of nanowave phononic 
devices including, e.g., mirrors, cavities, and monochro- 
matic sources. 

How long does it take for a particle to tunnel through 
a potential barrier? This is a question that has occupied 
physicists for decades and one for which there is still no 
definitive answer [3]. The Hartman effect (HE) states that 
the tunneling time becomes independent of the barrier 
length [4]. The independence of tunneling time on barrier 
length would imply arbitrarily large and indeed superlu- 
minal velocities for tunneling wave packets, if tunneling 
was in fact a propagation phenomenon. 

Phonon tunneling studies have also revealed phenom- 
ena related to the HE. Recent experiments on tunneling 

'Correspondence: arriaga@ifuap.buap.mx 

' Instituto de Fisica, Benemerita Universidad Autonoma de Puebla, Apartado 
Postal J-48, CP 7250 Puebla, Mexico 

Full list of author information is available at the end of the article 



acoustic waves have reported the breaking sound bar- 
rier [5,6]. Yang et al. found that inside a phononic band, 
the group velocity increases linearly with the sample 
thickness, a rather remarkable effect that is a signature 
of tunneling in quantum mechanics [6,7]. Villegas et 
al. have discussed the physical conditions under which 
the tunneling time for long-wavelength phonons through 
semiconductor heterostructures is independent on the 
system's size, i.e., the effect equivalent to the HE for elec- 
trons [8]. Experimental studies of HE in the context of 
the nanophononics have been carried out [9,10]. In these 
works, very short transit times in the stop bands have 
been measured, one acoustic equivalent of HE of electron 
tunneling through potential barriers. 

During the last decade, interest in achieving all-silicon- 
based opto- and microelectronics was highly stimulated 
by the discovery of the unique optical properties of porous 
silicon [11]. Porous silicon is known as a versatile material 
with applications in light emission, sensing, and photonic 
crystal devices. It is well known that the introduction of 
artificial spatial periodicity in the elastic properties of a 
system results in Brillouin zone folding. Such folding is 
often accompanied by the appearance of bandgaps in the 
phonon frequency spectrum [12,13]. In the last few years, 
this interest has been translated to porous silicon-based 
phononic systems [14-17]. Here, we study the tunneling 
times of acoustic phonon packets through a distributed 
Bragg reflector in porous silicon multilayer structures. 
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The paper is organized as follows: The 'Methods' 
section provides some fundamentals in both the long- 
wavelength model and the transfer matrix method. The 
main theoretical findings are presented in the 'Results 
and discussion' section. In particular in this section, from 
general principles such as invariance under time rever- 
sal, invariance under space reflection, and conservation of 
energy density flux, the equivalence of the tunneling times 
for both transmission and reflection is demonstrated. At 
the end of the paper, the main conclusions are given. 

Methods 

Long-wavelength model 

The one-dimensional energy density [18] is defined as 



u 



i 



pcof\u\ 2 - 



du* 
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, du 
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(1) 



where the first term in (1) represents the kinetic energy 
density, the second one represents the interaction energy 
density of the phonon field with itself, and the third one 
represents the strain energy density that accounts for 
the dispersive character of the oscillations. These terms 
depend on the atomic relative displacements u, the lin- 
ear mass density p, the phonon frequency at the center of 
the Brillouin zone &>r> the one-dimensional strain tensor 
du/dz, and the stress tensor a , which is equal to 



a = - P p 2 du/dz, 



(2) 



being /3 a parameter that accounts on the behavior of the 
bulk phonon dispersion relation. From Equations 1 and 2, 
we can obtain the one-dimensional equations of motion, 



and 
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The continuity equation for the energy density is 
defined as 



dH/dt + dj/dz = 0, 
where the energy density flux, /, is given by 
;' = -1/2 (adu*/dt + a*du/dt) . 



(5) 



(6) 



By substituting Equation 2 into Equation 6, we obtain the 
following expression for the energy density flux: 



; 



pP 2 (du* du du du* 
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(7) 



A very convenient method for calculating the propagation 
of harmonic waves through a system consisting of a finite 



number of uniform layers is the transfer matrix method 
[19,20]. The transfer matrix, which propagates the wave 
amplitudes across a j-th homogeneous layer is given by 



Mi 



cos (kjd,) sin (kjdj) 

a>Zj sin (kjdj) cos (kjdj) 



(8) 



where Zj = pjVj is the acoustic impedance, pi the mass 
density, vi the velocity of sound, and di the layer width. 

Results and discussion 

Equivalence of the transmission and reflection times in the 
tunneling of long-wavelength phonons 

In this paper, for the reason of complement, we give 
a demonstration of the equivalence of the transmission 
and reflection times (r t = T r ) in the tunneling of long- 
wavelength phonons [21]. Consider the general phonon 
scattering process as shown schematically in Figure 1. Let 
us consider that the phonon propagation is normal to 
the layer interfaces and adopt the continuum model valid 
for long-wavelength oscillations. In this section, we ana- 
lyze the properties of the transfer matrix in the tunneling 
process of long-wavelength phonons. 

invariance under time reversal 

Let us observe that if we take the complex conjugate of (3) 
and if t is replaced by —t, we get 



dt 2 



u*(z,-t) 



-a> 2 (z)u*(z, 



-t)-p( Z y 



dz 2 



u*(z, —t), 
(9) 




Figure 1 General stationary scattering configuration. The general 
stationary scattering configuration in one dimension. An arbitrary 
barrier is confined to the interval (—a, a). 
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provided only that f3(z) and a>r(z) are real functions. 
Observe that this equation has the same form as (3). 
Therefore, if u(z,t) is a solution of (3), then u*{z,—t) is 
also a solution. u*(z, — t) is often referred as the time- 
reversed solution. The behavior of the wave equation 
exhibited by (9) is called invariance under time reversal. 
For the stationary state, invariance under time reversal 
implies that if u(z) is a stationary-state wave function, 
then u* (z) is also one. 

The general solution of the wave equation for the system 
depicted in Figure 1 is given by 



ui(z) 



Ae ik ^ z + Be~ iklZ z < -a 

k 2 ) —a<z<a 
Ce ihz + De -ihz z> a 



(10) 



where An = ^ (co 2 Fi - to 2 ) /ft 2 (k 3 = J(a$ 3 - a?) / ft 

is the wave number in the left (right) part of the sys- 
tem. Using the transfer matrix, we can relate the wave 
amplitudes in the left and right parts according to 



where 
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On the other hand, the time-reversed solution of the wave 
equation is 



u 2 {z) 
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Hence, comparing Equations 10, 11, and 13, we get 

(5)-* 



where 
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The properties of the transfer matrix in the tunnel- 
ing of long-wavelength phonons can be obtained from 
Equations 11 and 14 



Af 
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Conservation of energy density flux 

To demonstrate the conservation of the energy density 
flux, we need to calculate this quantity in the left and right 
parts of the system. Using Equations 7 and 10, we obtain 
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where the superscript L (R) refers to the left (right) of 
the region (—a, a), Im denote imaginary part and Z\ = 
P\j3 2 k\ and Z 3 = p 3 p 2 k 3 . The conservation of the energy 
density flux is expressed as j L = ;' R , or equivalently 



M'lM* 



Zi 



I, 



(19) 



where M t and M* denote the transpose and conjugate 
matrix, respectively, and 



(o -0 



(20) 



Using Equation 19 and the condition (16), one further 
condition on matrix M can be added 



Z 3 

\M\ = MnM 22 - Mi 2 M 21 = 

Z\ 



(21) 



Invariance under space reflection 

Now, let us consider the invariance under space reflection 
of Equation 3. Because j3 2 {z), cd 2 {z), and io 2 {z) are even 
functions of z, another solution of the wave equation is 
obtained by replacing z by —z. We can immediately write 
the general solution of the wave equation as 

Ae~ ik i z + Be iklZ z> a 
u 3 (z) = { z,k 2 ) a > z > —a , (22) 

Q e -ihz + De ihz z < _ a 



and obtain the following relation 



M 



Substituting (23) into (14) gives us the condition 
MM = I, 



(23) 



(24) 



where / is the identity matrix. From this last equation, we 
obtain 

\M U \ 2 + \M l2 \ 2 = I and M n (M 12 + M* u ) = 0. 

(25) 

The transmission (reflection) amplitude Af (Af) for a 
phonon incident perpendicular to the barrier from right 
to left (the inverse process is described by D = 1, C = 
Af,A = 0,B = Af) and the transmission (reflection) 
amplitude A\ (A\) for a phonon incident from left to 



Lazcano etal. Nanoscale Research Letters 2014, 9:449 
http://www.nanoscalereslett.eom/content/9/1/449 



Page 4 of 7 



right (the direct process is described by A = 1,B = A\, 
C=A\,D = 0) are given by the following equations: 

A?(k 3 ,k 1 )=^ = \A?(k 3 ,k 1 )\e k *^\ 

A\ (* lf k 3 ) = -j- = |^ (h, h) I e<^\ (26) 



^ ( ^ fe) = = \ A L (kh h) | ^(faA), (27) 
Mn 

in terms of the modulus and phase. By applying conditions 
(16) and (21) to the previous relations, we obtain 



Af (k 3 , ft) = M = ^A\ (h, k 3 ) , 
Mn Zi 



A? (*3,*l) = -^Aj'(*i,A-3). 

M 2 i 



(28) 



(29) 



Finally, the transmission and reflection coefficients are 
given by 
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where t, r, and i denote transmitted, reflected, and inci- 
dent, respectively. From Equations 16, 27, 28, and 30, it is 
very easy to prove that T R (k 3 ) = T h {h) and "K R (fe) = 
lZ L (ki). From Equation 28, we obtain af(k 3 ,ki) = 
a\{ki,k 3 ), and from Equations 16, 26, 27, 28, and 29, we 
can deduce the relation 

P?(k 3 , ki) = ±it- ft{h, k 3 ) + 2a\{h, k 3 ). (31) 

If we consider that both materials in the left and right 
parts of Figure 1 are identical, i.e., (k\ = k 3 = k), 
then a t R (£ 3 ,£i) = (*Hh,k 3 ) = a t (k) and ^(k 3 ,h) = 
fi]:{k\,k 3 ) = P v (k). By using this relations, the phase is 



a t (k) = ±tt/2 + 6 r (k). 



(32) 



Finally, using Equation 32 is straightforward to prove that 
the transmission and reflection phase times given by 



dc/t 
dm ' 



and 

r da) ' 
are equal, 

T t = r r . 



(34) 



(35) 



(33) 



Falck et al. have obtained a similar result for symmetric 
scattering potentials [22]. 

Tunneling of acoustic phonons through a distributed 
Bragg reflector 

In this section, we apply the concepts discussed previ- 
ously to study a distributed Bragg reflector (DBR) formed 
by a finite number of periods, N, based on porous sili- 
con (PSi). The DBR is obtained by stacking periodically 
N times two constituent layers A and B with different 
porosity. The thicknesses of A and B layers are denoted 
by dA and dg, respectively. In PSi, the mass density is 
related with the porosity P via p = po (1 — P), being po 
the mass density of bulk silicon. The propagating veloc- 
ity of the longitudinal waves through PSi is related with 
the porosity as v = vq (1 — P) k , being vq the longitudinal 
velocity of the wave in bulk silicon, and A" is a param- 
eter. The acoustic impedance of layer A (B) is given by 
Za = Pava(Zb = PbVb)- The dependence on porosity of 
these two parameters, velocity and mass density, requires 
a very accurate control of the etching process. We grow 
the layers forming the DBR according to the methodol- 
ogy reported in [23]. The thicknesses of the PSi layers 
are dA = 2 ^m and dg = 1-65 pm, and their respec- 
tive porosities equal to Pa = 0.47 and Pg = 0.67. The 
parameters k = 0.56, po = 2,330 kg/m 3 , and vq = 
8, 440 m/s were determined by fitting the experimental 
results as is explained in [23]. In Figure 2, the theoret- 
ical (black solid line) and experimental (blue solid line) 
transmission coefficients (in dB), T = jt/ji, are plotted as 
a function of the acoustic phonon frequencies (in GHz). 
The theoretical transmission spectrum has been modeled 
using the transfer matrix method and includes the effect 
of the sample, transducers, and liquid coupling the trans- 
ducers to the sample. All these components are included 
writing their corresponding transfer matrix like appearing 
in (8) but using their respective parameters. The exper- 
imental transmission has been measured using a vector 
network analyzer (VNA) using two piezoelectric trans- 
ducers according to the experimental setup described in 
[23]. We can observe that in the range of frequencies 
reported here, the DBR shows two acoustic gaps, corre- 
sponding to the first and second acoustic gap, and cen- 
tered around = 0.71 and = 1.42 GHz, respectively. In the 
inset, we show the calculated phonon transmittance for 
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Figure 2 Transmission coefficients. Theoretical and experimental transmission coefficients (in dB) plotted as a function of the acoustic phonon 
frequencie (in GHz) for the DBR, with N = 18. The inset shows some details of the theoretical transmission coefficient around the first acoustic stop 
band. 



frequencies around the zone center of the first minigap. 
In particular, note that at frequencies fi = 0.54 GHz and 
fg = 0.88 GHz, corresponding to the edges of the gap, 
T = 0 dB. We can observe a noticeable difference in 
the fundamental bandgap between the experimental and 
the theoretical spectrum. The experimental bandgap has a 
depth of approximately 50 dB which is less than the mod- 
eled value of approximately 200 dB. However, this is due 
to the experimental limitation on the setup used, and it is 
only attributed to the noise of our VNA. 

The behavior of the vibrational energy density % (1) as a 
function of the normalized distance, along the axis of the 
DBR, is shown in Figure 3a. We observe that the vibra- 
tional energy density is a one piecewise constant function. 
Let us now investigate the behavior of H. at the reso- 
nances. At the frequencies ft = 0.54 GHz (blue solid line) 
and /r = 0.88 GHz (red solid line), the vibrational energy 
density exhibits the maximum value in the center region 
of the DBR. However, these modes are characterized by 
complete transmission of the vibrational energy through 
the system. The one-dimensional stress tensor a and the 
relative displacements u are shown in Figure 3b,c, respec- 
tively. In these figures, we have included a schematic rep- 
resentation of the structure, representing with dark (light) 
regions the layers with low (high) porosity. We observed 
from Figure 3c that the relative displacement for the mode 
appearing at /l = 0.54 GHz (/r = 0.88 GHz) is localized 
in the low (high) porosity regions, which is consistent with 
the agreement used to refer at the valence (dielectric) and 
conduction (air) band used in semiconductors (photonic 



crystals). From these figures, it is clearly observed that the 
phonon modes propagate through the structure. 

The phonon modes with energies within the first mini- 
gap appear marked by the black solid line in Figure 3b,c. 
We observe that at the frequency fc = 0.70 GHz, the 
amplitude of the wave shows an abrupt decay along the 
axis of the DBR. This qualitative behavior is very similar 
for the other physical quantities. 

Tunneling times 

The dwell time, t\), for acoustic phonons inside the DBR is 
calculated by the following equation [21]: 
i_i 

to = ~. (36) 
h 

where H is the total energy in the interval (0,L), being L 
the total length of the DBR. This energy is obtained by 
integrating Equation 1, i.e., 

H= f Hdz. (37) 
Jo 

The free time, t[, which is the time associated to the transit 
of the sound pulse along a distance L with velocity, Vf, is 
simply defined by 

t f = (38) 

Vf 

being V{ the phase velocity. 

In Figure 4, the tunneling times as a function of the 
phonon frequency (in GHz) are depicted for the acous- 
tic phonons around the second acoustic gap. In Figure 4a, 



Lazcano etal. Nanoscale Research Letters 2014, 9:449 
http://www.nanoscalereslett.eom/content/9/1/449 



Page 6 of 7 



(a. u.) 




f c =0.70GHz 

J J~ n f L =0.54GHz- 

I— T~| - f R =0.88GHz_ 


* 5 






0 




_ 1.2 






«. 0.6 






"5" 

0.0 







3 



iiiiiiii iiiiii i 




Figure 3 The vibrational energy density, the one-dimensional stress tensor and the relative displacements, (a) The vibrational energy 
density H (in a.u.), (b) the one-dimensional stress tensor a (in a.u.), (c) and the relative displacements u (in a.u.) in the structure, plotted as a function 
of the normalized distance along the axis of the DBR. The distance is normalized by the factor Ic = + dg. 



we can observe that the dwell time (36) increases near 
the gap. The behavior at such frequencies is the same as 
if the phonon was trapped for a long time in the spatial 
region occupied for the structure before being transmit- 
ted. In Figure 4b, we observe that, in the stop bands, 



the transmission time (33) is shorter than free time (38), 
though a large noise is observed due to the coexistence of 
the liquid modes. The very short transit times in the stop 
bands are the acoustic equivalent of the Hartman effect 
[4] of electrons tunneling through potential barriers [8,9]. 
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Figure 4 Tunneling times, (a) The dwell time to (in u.s) and (b) the transmission time r t and free time ff (in its) plotted as a function of the acoustic 
phonon frequencies (in GHz), (c) Theoretical and experimental transmission coefficients (in dB) plotted as a function of the acoustic phonon 
frequencies (in GHz). 
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Conclusions 

In this paper, we studied tunneling times of acoustic 
phonon packets through a distributed Bragg reflector 
made of porous silicon layers. Under the assumption that 
the long-wavelength approximation is valid, and from 
general principles of symmetry and conservation, we 
report an explicit demonstration of the equivalence of 
the transmission and reflection times in the tunneling of 
long-wavelength phonons. Calculations of the vibrational 
energy density and the vibrational energy stored within 
the structure allows a better visualization of the physi- 
cal phenomena occurring in this system. The description 
of the stress and strain fields complements the energetic 
description. We report the possibility that a phenomenon 
called Hartman effect appears in porous silicon multilayer 
structures, an acoustic equivalent of Hartman effect of 
electrons tunneling through potential barriers. The results 
of this study could be useful for the design of acoustic 
devices. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

All authors contributed to the analytical calculations. ZL performed the 
experimental measurements and fabrication of the porous silicon DBR. All 
authors contributed to the writing of the manuscript. All authors read and 
approved the final manuscript. 

Acknowledgements 

DV and RPA acknowledge the hospitality of Benemerita Universidad Autonoma 
de Puebla, Mexico. DV acknowledges the support of S.R.E. (Mexico). JA 
acknowledges Conacyt, Mexico for the partial support under grant no. 1 67939. 

Author details 

1 1nstituto de Fisica, Benemerita Universidad Autonoma de Puebla, Apartado 
Postal J-48, CP 7250 Puebla, Mexico. 2 Departamento de Fisica, Universidad 
Central 'Marta Abreu' de Las Villas, CP 54830 Santa Clara, Cuba. 3 Universidad 
Autonoma del Estado de Morelos, Ave. Universidad 1 001 , CP 62209 
Cuemavaca, Mexico. 



8. Villegas D, de Leon-Perez F, Perez-Alvarez R: Tunneling time for 
phonons: dependence on the system's size. Phys Status Solidi (b) 2005, 

242:1767-1770. 

9. Huynh A, Lanzillotti-Kimura N, Jusserand B, Perrin B, Fainstein A, 
Pascual-Winter M, Peronne E, Lemaitre A: Subterahertz phonon 
dynamics in acoustic nanocavities. Phys Rev Lett 2006, 97:1 15502. 

1 0. Lanzillotti-Kimura ND: Dinamica de fonones acusticos en multicapas 
nanometricas. PhD thesis Institute Balseiro, SC Bariloche, Argentina 2009. 

1 1 . Ghulinyan M, Otton CJ, Gaburro Z, Bettotti P, Pavesi L: Porous silicon 
free-standing coupled microcavities. Appl Phs Lett 2003, 82:1 550-1 552 

1 2. Parsons LC, Andrews GT: Observation of hypersonic phononic crystal 
effects in porous silicon superlattices. Appl Phys Lett 2009, 95:241 909. 

1 3. Aliev GN, Goller B, Kovalev D, Snow PA: Hypersonic acoustic mirrors 
and microcavities in porous silicon. Appl Phys Lett 2010, 96:124101. 

1 4. Da Fonseca PJM, Saurel JM, Foucaran A, Massone E, Talierco E, Camassel J: 
Acoustic microscopy investigation of porous silicon. Thin Solid Films 
1995, 225:155. 

1 5. Andrews GT, Zuk J, Kiefte H, Clouter MJ, Nossarzewska-Orlowska E: Elastic 
characterization of a supported porous silicon layer by Brillouin 
scattering. Appl Phys Lett 1 996, 69:1 2 1 7. 

1 6. Thomas L, Aliev GN, Snow PA: Hypersonic rugate filters based on 
porous silicon. Appl Phys Lett 201 0, 97:1 73503. 

1 7. Aliev GN, Goller B, Snow PA: Porous silicon bulk acoustic wave 
resonator with integrated transducer. Nanoscale Res Lett 201 2, 7:378. 

1 8. de Leon-Perez F, Perez-Alvarez R: Phonon propagation in nonpolar 
semiconductor heterostructures. Phys RevB 2001, 63:245304. 

1 9. Kavokin A, Malpuech G, Di Carlo A, Lugli P: Photonic Bloch oscillations 
in laterally confined Bragg mirrors. Phys Rev B 2000, 61 :441 3 

20. Karabutov AA, Kosevichc Y, Sapozhnikov OA: Bloch oscillations of an 
acoustic field in a layered structure. Acoust Phys 2013, 59:158-169. 

21 . Villegas D, de Leon-Perez F, Perez-Alvarez R: Tunneling time of 
long-wavelength phonons through semiconductor 
heterostructures. Phys Rev B 2005, 71:035322. 

22. Hauge EH, Falck JP, Fjeldly TA: Transmission and ref letion times for 
scattering of wave packets off tunneling barriers. Phys RevB 1987, 
36:4203. 

23. Lazcano Z, Arriaga J, Aliev GN: Experimental and theoretical 
demonstration of acoustic Bloch oscillations in porous silicon 
structures. J Appl Phys 20 1 4, 11 5: 1 54505. 

f N 

doi:1 0.1 186/1 556-276X-9-449 

Cite this article as: Lazcano et ah Tunneling times of acoustic phonon 
packets through a distributed Bragg reflector. Nanoscale Research Letters 
2014 9:449. 



Received: 1 0 May 201 4 Accepted: 1 2 August 201 4 
Published: 29 August 2014 



References 



Balandin AA: Nanophononics: phonon engineering in nanostructures 
and nanodevices. J Nanosci Nanotechnol 2005, 5:1015-1022. 
Narayanamurti V, Stormer HL, Chin MA, Gossard AC, Wiegmann W: 
Selective transmission of high-frequency phonons by a superlattice: 
the "dielectric" phonon filter. Phys Rev Lett 1 979, 43:201 2-201 6. 
Winful HG: Delay time and the Hartman effect in quantum tunneling. 
Phys Rev Lett 2003, 91 :260401 . 

Hartman TE: Tunneling of a wave packet. J of Appl Phys 1 962, 

33:3427-3433. 

Robertson W, Ash J, McGaugh J: Breaking the sound barrier: tunneling 
of acoustic waves through the forbidden transmission region of a 
one-dimensional acoustic band gap array. Am J of Phys 2002, 70:689. 
Page J, Suxia Y, Cowan M, Zhengyou L, Chan C, Ping S: 3D Phononic 
cristals. In Wave Scattering in Complex Media: from Theory to Applications. 
Edited by van Tiggelen B, Skipetrov S. Amsterdam: Kluwer Academic; 
2003:283-307. 

Yang S, Page J, Liu Z, Cowan M, Chan CT, Sheng P: Ultrasound tunneling 
through 3D phononic crystals. Phys Rev Lett 2002, 88:1 04301 



Submit your manuscript to a SpringerOpen 
journal and benefit from: 

► Convenient online submission 

► Rigorous peer review 

► Immediate publication on acceptance 

► Open access: articles freely available online 

► High visibility within the field 

► Retaining the copyright to your article 

Submit your next manuscript at ► springeropen.com 



